Is J enough? Comparison of gravitational waves emitted along the total angular 
momentum direction with other preferred orientations 
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The gravitational wave signature emitted from a merging binary depends on the orientation of 
an observer relative to the binary. Previous studies suggest that emission along the total initial 
or total final angular momenta leads to both the strongest and simplest signal from a precessing 
compact binary. In this paper we describe a concrete counterexample: a binary with mi/m2 = 4, 
ai = 0.6x = —a2, placed in orbit in the x, y plane. We extract the gravitational wave emission along 
several proposed emission directions, including the initial (Newtonian) orbital angular momentum; 
the final (~ initial) total angular momentum; and the dominant principal axis of (L(a-£/6))^^. Using 
several diagnostics, we show that the suggested preferred directions are not representative. For 
example, only for a handful of other directions (p < 15%) will the gravitational wave signal have 
comparable shape to the one extracted along each of these fiducial directions, as measured by a 
generalized overlap (> 0.95). We conclude that the information available in just one direction (or 
mode) does not adequately encode the complexity of orientation-dependent emission for even short 
signals from merging black hole binaries. Future investigations of precessing, unequal-mass binaries 
should carefully explore and model their orientation-dependent emission. 
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I. INTRODUCTION 



Numerical simulations of merging compact binaries 
produce a gravitational wave signature that depends 
on the orientation of the binary relative to the line of 
sight. For nonprecessing comparable-mass binaries, this 
orientation dependence can be wcU-modclcd as purely 
quadrupolar emission from conjugate and orthogonal 
I = \m\ = 2 modes. In this approximation, only one 
complex-valued function is needed to characterize the 
gravitational wave signature from a binary: either the 
{l,m) — (2,2) mode itself or the strain extracted along 
the binary angular momentum axis. This approach has 
been broadly adopted when comparing nonspinning nu- 
merical relativity simulations to one another [l| and to 
post-Newtonian and other Q approximations; when 
building hybrid waveforms that join systematic post- 
Newtonian approximations to numerical relativity 0, Q ; 
when constructing phenomcnological approximations to 
numerical relativity waveforms [9nlh]: and when search- 
ing for the gravitational wave signature of merging com- 
pact binaries with interferometric detectors [l2|. For 
similar reasons, the gravitational wave signature from a 
nonprecessing unequal-mass binary can also be approxi- 
mated by h^^x{t, J), the radiation along the total angular 
momentum axis J. 

Motivated by this approximation, some have proposed 
that radiation along "the total angular momentum di- 
rection" is a fiducial example of the gravitational wave 
signature from a binary [l3|. In general, the total angu- 
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lar momentum direction j is not conserved during the 
inspiral and coalescence [isl - ITTj . To a good approx- 
imation, however, the orbital angular momentum often 
precesses around the total angular momentum . Dur- 
ing the inspiral, the binary radiates predominantly along 
the instantaneous orbital angular momentum [ij; fl8l - [20| , 
which precesses around the total angular momentum. In 
this case, previous direct and indirect calculations sug- 
gest that that strong, simple emission will occur along 
the average direction, which for simple precession corre- 
sponds to the total angular momentum [l^, [2lj . In this 
paper, we provide an explicit counterexample, for which 
J neither corresponds to nor produces radiation simi- 
lar to the peak emission direction. We contrast J with 
the time-dependent and signal-averaged directions from 
(^L(^aLb)} 2l|, which better reflects the peak emission di- 
rection. Finally, we explicitly demonstrate that higher 
harmonics and precession prevent a significant fraction 
of orientations from being well-fit as the linear super- 
position of only two basis waveforms in an "antipodal 
decomposition." The orientation-dependent emission of 
precessing waveforms is unavoidably more complicated 
than the response of their nonprecessing counterparts. 



A. Context 

Merging binaries imprint their properties on the grav- 
itational waves radiated to infinity in all directions 
In practice, however, gravitational wave detectors have 
access to only one line of sight, sometimes only one po- 
larization. Under the right circumstances, two physically 
distinct configurations can produce very similar emission 
in the detectors' sensitive band. These near- deg eneracies 
complicate source parameter estimation [22h27{ and thus 
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the interpretation of gravitational wave measurements for 
astrophysics. 

Conversely, a single source can produce highly dissim- 
ilar signals along different lines of sight. By contrast, the 
signal template used in real searches for nonprecessing 
binaries in gravitational wave data assumes a universal 
form for all orientations: a linear superposition of two 
basis functions, the "cosine" and "sine" chirps [l2|. For 
nonprecessing binaries, this signal model is physically 
well-motivated [see the Appendix]; the relative weights 
of the modes are in one-to-one relation to the source in- 
clination relative to the line of sight. In this paper we 
use a concrete example where this decomposition fails. 
For the example described in this paper, more than two 
functional degrees of freedom are needed to reproduce 
the gravitational wave signal accurately for all lines of 
sight. 



II. ORIENTATION-DEPENDENT EMISSION 

We simulated the evolution of a binary black hole 
(BBH), with total mass M = mi + m2 = 1.0 and mass 
ratio toi/to2 = 4, initiated on the a;-axis, with initial 
separation d = 9 and initial dimensionless spin vectors 
ai = 0.6i = — a2. For reference, the initial and final 
total angular momenta are 



Ji ~ 0.36a; 
Jf = 0.27x 



0.602 
0.036y 



0.49z 



(1) 
(2) 



Initial data was evolved with Maya, which was used in 
previous BBH studies [28l - l36j . The grid structure for 
each run consisted of 10 levels of refinement provided 
by CARPET [33, a mesh refinement package for CACTUS 
[38| . Sixth-order spatial finite differencing was used with 
the BSSN equations implemented with Kranc [s^. The 
outer boundaries are located at 307M. Each simulation 
was performed with a resolution of M/140 on the finest 
refinement level, with each successive level's resolution 
decreased by a factor of 2. In the text, we present re- 
sults for the fiducial r = 90 extraction radius and for 
this resolution. The comparatively large differences we 
describe exist at all large extraction radii and simulation 
resolutions we have explored. Similarly, though we adopt 
fiducial snapshot parameters (time and mass) to generate 
the contours in Figures [1] and [2l we find similar behavior 
at neighboring masses and times. 

This simulation has markedly anisotropic emission. 
Figure [T] shows contours of 1-041 (solid) and dE/dtdft oc 
J ipidtl'^ (dashed), superimposed with black dots indi- 
cating the total angular momentum direction and a blue 
dot indicating the dominant eigenvector of (^Lf^aLb)} j2]| . 
This tensor is defined by 

/r r \ _ Co.dfJdnr,{f)haLbMf)/Sh 
J-oodf J dVL\^i{f)\^ / Sh 



where Sh is the noise power spectrum of a candidate in- 
terferometer. None of these three directions corresponds 
to the directions of peak transverse tidal force or radiated 
power. 

For clarity we only show contours for a single time. 
Over the inspiral, the beampattern evolves slightly, with 
the orientation of strongest -04 and dE / dt roughly follow- 
ing the path of our preferred orientation (blue line) . For 
all times simulated, the largest -04 and dE/dtdfl is never 
along the total angular momentum. 

For comparison, in the bottom panel of Figure [1] we 
show contours for a simulation with comparable mass ra- 
tio in the absence of spin. Though vertically symmetric, 
the functional form of "04 and dE / dtdU, is largest at spe- 
cific angle that evolves along with the (retarded) sepa- 
ration vector. For this particular timeslice, the peak oc- 
curs at (6*, 0) ~ (0.5, 1.5), significantly offset from the an- 
gular momentum axis. This shape rotates at the orbital 
period through merger. By contrast, all three preferred 
orientations {L, J and V , the principal eigendirection of 
l^L^^^L^'^) always lie along the z axis. 

For data analysis purposes, a more relevant 
orientation-dependent quantity is the coherent signal-to- 
noise accessible to a particular detector network along 
that line of sight. Under the simplifying assumptions 
that two identical interferometers misaligned by 7r/4 
intercept gravitational wave emission perpendicular to 
the line of sight, an idealized orientation-dependent 
signal-to-noise ratio p{n) is given by [2lll40l|^ 



df 1^4 (/,n) I 
Sh (27r/)4 



(4) 



where Sh is the noise power spectrum of each of the two 
interferometers. (This expression is simply the sum of 
the signal to noise due to each detector {p\ + P2), using 
"04 = d^{h+ — ihy^) to relate the sum to a quantity nat- 
urally produced by simulations.) The dotted contours in 
Figure [5] show contours of constant p, adopting a sim- 
ple analytic model for the initial LIGO noise curve [4l| 
and a fiducial total mass M = lOOM© for the binary 
to scale our (dimenionless) simulations. Again, neither 
the total nor (Newtonian) orbital angular momentum di- 
rection seems significant with respect to the detectable 
signal power. 

The contours of and its peak location inevitably are 
qualitatively similar to contours of dE/dfdfl, the energy 
radiated per unit frequency. The expression for p^(M) 
can be rewritten as a suitable frequency-weighted aver- 
age of dE/dfdil cx I -04 //p. For a suitable choice of /, 
contours in the neighborhood of the maximum are qual- 
itatively similar. 



^ For comparison, the orientation-averaged signal to noise p de- 
scribed in O'Shaughnessy et al. [4(1 is precisely the orientation 
average (p^) of our signal- to- noise ratio over all emission direc- 
tions. 
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Two- or more detector networks are generally sensi- 
tive to gravitational wave polarization. Our cartoon two- 
detector network has equal sensitivity to both polariza- 
tions. Such a network can automatically project out left- 
or right-handed signal. Motivated by sensitivity to helic- 
ity, we define two signal to noise ratios p± that measure 
sensitivity to right and left handed cmisison, respectively: 



pI 



p'- 



df \Mf,n)\' 

df \Mf,n)\' 
_^S„ (27r/)4 



(5) 
(6) 



At each instant, right-handed emission occurs along the 
orbital angular momentum direction, left-handed op- 
posite. These expressions encode the contributions of 
right-handed and left-handed emission to the total SNR 
= p'^_ + p\. The solid curves in Figure [2] are con- 
tours of p+. Their shape and similarity to contours of 
p demonstrate that each of the two local extrema of p^ 
is produced near-exclusively by emission with a single 
handedness. Neither local extrema corresponds to the 
direction of J or the (Newtonian) orbital angular mo- 
mentum. By contrast, the nearly-antipodal extrema do 
correspond to the orientation of the principal eigenvalue 
of (L(aL(,))^^ dll, shown as a red dot in Figure^] 

To summarize, for our fiducial binary, we have exam- 
ined the gravitational wave signature along three com- 
mon choices for the preferred orientation: the princi- 
pal axis of (L(a^b))jv/ as well as the (Newtonian) or- 
bital and total angular momentum directions. Neither 
angular momentum direction corresponds to a direction 
extremizing strain, energy flux, or detectable signal to 
noise. By contrast, the direction identified by (L(aib))^ 
nearly corresponds to a local maximum of the signal-to- 
noise p^ . We anticipate and will show below that the 
direction identified by (^Lf^aLb)) j^j is "more representa- 
tive" of the orientation-dependent emission, in that a 
greater fraction of the orientation-dependent emission 
can be well-fit by it. However, none of the three direc- 
tions can adequately fit all the complexities of this bi- 
nary's orientation-dependent emission, even for one rep- 
resentative mass Af = IOOMq. 



III. CONTRASTING DIFFERENT 
ORIENTATIONS 

Most variations in the gravitational wave signature or 
amplitude along different emission directions are too sub- 
tle to be distinguished by gravitational wave networks at 
the amplitude of astrophysically plausible signals. Gener- 
ally, distinguishability of signals depends strongly on the 
detector network, signal model, and source location rel- 
ative to the earth. Rather than address this complicated 
issue in general, we adopt a simple physically-motivated 
proof of concept calculation. The complications associ- 
ated with real detector networks will, if anything, further 
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FIG. 1: Time-domain contours of 1 1/)4 1 (solid) anA dE / dtdO. 
(dashed), evaluated at the peak of the {l,m) — (2,2) mode. 
Both quantities are normalized to 1 at S = = 0. Solid 
black dots indicate the location of the initial and final total 
angular momenta. The solid blue dot indicates the preferred 
axis identified by ^L(t,Lj))^ at this time. The dark blue path 
indicates how this preferred orientation evolves with time, 
several hundred M before and 30Af after this reference time. 
The top panel corresponds to the fiducial spinning binary. 
The bottom panel shows a binary with comparable mass ratio 
(mi/m2 = 4) but without any spin. 
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break any symmetries connecting the natural emission 
direction to J. 

First and foremost, we introduce a complex overlap 
and its natural maximization in time and polarization, 
a tool which efficiently eliminates two extrinsic param- 
eters of generic signals even in the presence of higher 
harmonics. Though we adopt this complex overlap for 
technical convenience and its explicit relationship with 
(^(aib)), all results described below also follow from the 
standard inner product on a single polarization at a time 
[Eq. ([TT]) : see ^42, 43]]. Technically, subtle differences 
arise due to the number of data streams and the order of 
maximization over extrinsic parameters: time and polar- 
ization (two-detector complex) or time and phase (single- 
detector match). Practically speaking, however, the two 
expressions often return comparable results, particularly 
when one of the two signals is nearly circularly polar- 
ized as will be the case below. As an example, in Fig- 
ure [3] we compare the complex match P to the "typical" 
and "minimax" matches [42|. In regions of high com- 
plex overlap, the two are qualitatively similar. In regions 
with low complex overlap, the "typical" and "minimax" 
matches differ substantially, indicating the two linear po- 
larizations along that line of sight are distinguishable and 
comparable in magnitude. 



A. Complex overlap 

For each emission direction h, we once again assume 
two identical interferometers are placed, oriented by 7r/4 
relative to one another and perpendicular to the incident 
signal. We coherently compare the (noise-free) signal ex- 
pected along any pair of orientations with a complex 
inner product motivated by the detector's noise power 
spectrum pTj . 

For our purposes, numerical relativity simulations take 
as inputs binary black hole parameters and desired line 
of sight (denoted by A) and return the Weyl scalar V'4(i), 
a complex-valued function of time evaluated along that 
line of sight. For any pair of simulations and lines of 
sight, we compare '04 and ip4 by a complex overlap 



FIG. 2: Detector-weighted contours of p (dotted) and p+ 
(solid), evaluated at M — lOOi\/0, plotted versus the simu- 
lation's asymptotic spherical polar coordinates (i.e., aligned 
with the initial orbital angular momentum). Solid black dots 
indicate the location of the initial and final total angular mo- 
menta. The solid red dot indicates the preferred axis identi- 
fied by {^L(aLb)) _^ [Eq. (jAl|) ] at this time; nearly identical 
preferred axes are identifies by (Lj^L;,))^^ and (L(;,L[,))^^_. 
As in Fig. [TJ the solid blue path indicates the evolution of 
(L(t,Lj)) versus time. The dashed green path indicates the 
orientation of the (Newtonian, coordinate) orbital angular 
momentum oc r x v. The top panel corresponds to the fiducial 
spinning binary, while the bottom panel shows an nonspinning 
binary with comparable mass ratio (g = 4). In the bottom 
panel, the preferred orientation identified by {L{aLb)) ^ is 
slightly above the 6 = line; the two directions identified by 
(^Li^aLh)) are not precisely antipodal. 



P(A,A') ^ 
(A,B) ^ 



\r'il)i\\r-4}i\ 



df 



(7a) 

i(/)*B(/) (7b) 



where Sh is a detector strain noise power spectrum. In 
this and subsequent expressions we used unprimed and 
primed variables to distinguish between the two wave- 
forms being compared, involving potentially distinct pa- 
rameters A' and lines of sight n' . For simplicity and to 
avoid ambiguity, in this paper, we adopt a semianalytic 
model for the initial LIGO sensitivity [4l|. As with the 
single-detector overlap, the complex overlap can be max- 
imized over the event time and polarization {t^^'c) and 
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by a simple fourier transform: 



Pr, 



maxt^^^jP 
1 

l'/'4||V'4'| 



(8) 



df 



MfTMf) 



For a network with isotropic sensitivity to both gravi- 
tational wave polarizations,^ this complex-valued expres- 
sion naturally characterizes the ability of the network 
to differentiate signals. The real part of this expression 
corresponds to a coherent sum of the conventional over- 
laps of the two gravitational wave polarizations generated 
from ^p•. 



{A\B) 



{h+\h'+) + {h^\h'^) 
-iif)/{2nff 

2 / dfA*B/Sh 



(9) 
(10) 

(11) 



where {A\B) is the standard single-detector overlap of 
gravitational wave strain. This expression also relies on 
the gauge-invariant output of gravitational wave simula- 
tions, without any complications associated with conver- 
sion ?/;4 — )■ h{t). Critically, this expression naturally re- 
spects rotational symmetry around the line of sight: if the 
two simulations are rotated by 7,7', then 1^4^ — >■ e~^*'^'!/'4 
and similarly, so P changes by a pure phase: 

P{RX,R'X') = e-2^('^'-'^'P(A, A') 

Equivalently, the absolute value |P| corresponds to choos- 
ing the (physically arbitrary) orientations of the two sim- 
ulations such that P is real. Similarly, the complex over- 
lap interacts transparently with a decomposition of the 
NR signal into (complex) spin-weighted spherical har- 
monics. For example, as with the conventional single- 
detector signal, a network with equal sensitivities to two 
polarizations has a network Fisher matrix along the line 
of sight proportional to 



Fab = (9aV'4, dblpA:) 



(12) 



For rotations, the Fisher matrix can be reorganized as an 
action of SU(2) generators on 1^4 [2l| : 



(13) 



This action can be evaluated algebraically using the 
known properties of SU(2). 

Qualitatively speaking, the complex overlap is con- 
siderably more discriminating than a single-polarization 



^ Currently-planned gravitational-wave detector networks arc pri- 
marily sensitive to one polarization for most sky locations 
l44l I45II . Future gravitational wave networks will try to achieve 
comparable sensitivity in both pol arizations, to better measure 
source inclination and distance |44||. 
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FIG. 3: Comparison of overlap for complex versus 
standard match: These curves compare three different cri- 
teria to distinguish between emission along z and n{9,<j)), 
shown for 6 = 40.5 deg (black) and 53deg (blue). The solid 
curve is the normalized complex overlap \P{z,n)\. The dot- 
ted curve and dashed curves are the "typical" and "minimax" 
match constructed from h+,hx along the two lines of sight 
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overlap. As an example, because the complex overlap 
naturally distinguishes between helicities, the complex 
overlap between the z and —z directions of an inspiralling 
nonprcccssing binary is nearly 0. By contrast, a single- 
detector overlap would be nearly unity, as one polariza- 
tion cannot distinguish between left- and right-handed 
signals. Finally, in the complex overlap, maximization 
over the physically irrelevant polarization angle 7 corre- 
sponds to an absolute value. 



B. Area of high overlap 

For each of the three reference orientations identified 
above - ii = z, the initial (Newtonian) orbital angular 
momentum direction; J, the nearly-conserved final an- 
gular momentum direction; and V, the preferred orien- 
tation identified by l^L^^aLb)) j^jj^ - we compute Pmax{v, n) 
everywhere on the sphere. 

For comparison, as indicated by the red curve in Fig- 
ure m for a nonprecessing binary the emission along 
+J = z ^ V alone can fit about one third of all emis- 



^ For a single-detector overlap, maximization over polarization re- 
quires diagonalizing a quadratic form for basis waveforms l42l . 
Physically, this process corresponds to choosing a detector orien- 
tation aligned with the natural principal axes of the waveform, 
projected on the plane of the sky. When multiple harmonics 
are present, maximization of single-detector overlaps becomes 
an arduous process; see, e.g., the Appendix of McWilliams et al. 
[4^ . Finally, for the single-detector overlap, efficient maximiza- 
tion over polarization is generally incompatible with fast (fou rier- 
transform- based) determination of the coalescence time |42|| . 
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FIG. 4: Area of high overlap: Fraction p of tiie spiiere 
for wliicii tiie overlap Pmax{v,fi) is above a given threshold 
P, versus the threshold, for v = z (blue), J (green), and the 
preferred axis of (^L(^a.Lf,)) (black). For comparison, the thin 
red curve shows the corresponding result for a binary radiat- 
ing only conjugate and orthogonal (2, 2) and (2, —2) modes 
and a reference axis v = z; see the appendix for a derivation. 
For this idealized nonprecessing binary, more than a third of 
all orientations have better than P = 0.95 (complex) match 
with the (2, 2) mode; all are well-fit with an antipodal decom- 
position. For the three reference orientations used here, the 
gravitational wave signal in that direction closely resembles 
the gravitational wave signal for only a small neighborhood 
surrounding that point. Results for simulations are estimated 
using ~ 10'' candidate viewing directions, each of which in- 
cludes alH < 4 contributions to ?/)4(?i). 



sion orientations, weighted equally in area; and about 
48%, weighted uniformly in (i.e., by equal probability 
of detection in the local universe). Summing over both 
antipodal directions, almost all detected emission from a 
nonprecessing binary is well- fit by either of 'ip4{±J). 

By contrast. Figure |4] shows that each of these direc- 
tions provides a good fit (P > 0.95) only for a small 
area surrounding each reference point. The gravitational 
wave signal along the (Newtonian) orbital and total an- 
gular momentum directions (i, J) are particularly non- 
representative, reproducing signals emitted only in im- 
mediately neighboring directions. 



C. Failure of antipodal reconstruction 

Existing gravitational wave searches for nonprecessing 
binaries, however, do not merely fit each data set against 
two candidate signals per simulation. Motivated by the 
simple functional form of mass-quadrupole-dominated 
emission, real searches project the data onto the two- 
dimensional subspace spanned by those two candidates."^ 



When the candidate signals are associated with two dif- 
ferent antipodal emission directions (±"0) we will denote 
this procedure an antipodal decomposition. 

Algebraically, an antipodal decomposition estimates 
the Weyl scalar along any direction via a projection to 
the subspace spanned by the two antipodes: 



\N) 
\S) 
Pn 

PS 
iciest 



M-i)/\M-v)\ 

(5, ip4{n))_ 
[{pN - {S, N)ps) \N) + {ps - {N, S)pn) 



(14) 
(15) 

(16) 



P 



\pn\^ - 



1-\{N,S)\' 
(7A4, (1 - V)^4) 



\Ps\ 



2Re(iV, S)pNP*s 



1-\{N,S)\ 



(18) 



Practically, an antipodal projection of this kind is useful 
if the relative loss of SNR is small (i.e., Prem/p < 0.05). 

For nonprecessing signals, the success of an antipodal 
approximation along nearly any axis is guaranteed by pa- 
rameter counting. Comparable-mass nonprecessing bina- 
ries are dominated by just two modes (the {I, m) = (2, 2) 
and (2, —2) modes), with two coefficient functions '042±2- 
As any antipodal basis arc in their span, an antipodal de- 
composition just corresponds to a change of basis. For 
nonprecessing binaries, however, the choice v = J pro- 
duces added benefits: the antipodes are local extrema of 
?/'4|, dE/dtdO,, and p^ . In this frame, trivially the angu- 
lar dependence of pN,s{n) decomposes into roughly two 
spin- weighted harmonics, to the extent that higher har- 
monics can be neglected. Finally, as noted previously this 
antipodal decomposition is physical: the ratio Pn/ Ps is 
in one-to-one relation with the inclination. 

By constrast, for precessing binaries the signal space 
generally has at least 5(= 2x2-1-1) basis vectors and thus 
has 5 complex (10 real) dimensions [13]. The precession 
of instantaneous (2, 2) and (2, —2) modes aligned with 
L projects into other I = 2 modes in the global frame. 
As a result, the "remaining signal power" prem left af- 
ter an antipodal decomposition is generally nonzero and 
often significant, if either precession occurs or if higher 
harmonics contribute significantly to overall radiation. 

Nonprecessing binaries also emit a significant fraction 
of their total (band-limited) signal p^ in higher harmon- 
ics, at high mass HO] or mass ratio 0, 13]. However, 
comparable-mass nonprecessing binaries are often still 
well-approximated as antipodal: for most masses, they 
emit most strongly in two complex-conjugate I = \m\ 



* Their real and imaginary parts arc known as the "cosine" and 
"sine" chirps. For sufficiently low mass, these two modes are 



nearly orthogonal. Their orthogonality can be exploited to effi- 
ciently maximize in time and "coalescence phase" [T3l . At high 
mass too few cycles exist to guarantee orthogonality of the (2, 2) 
and (2, —2) modes in general. 
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0.39 


0.44 
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0.01 
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0.009 
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0.045 
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-3 


0.037 


0.045 
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4 


0.005 


0.007 
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0.005 


0.007 



TABLE I: Multiple harmonics contribute significantly 

to the orientation-averaged signal power, for our reference 
simulation and mass (IOOM0). This table refers to mode 
amplitudes in a frame aligned with the prinicipal axis of 
(^L(^aLb)) j^j- The second and third columns show results for 
our two fiducial simulations, both with q — 4 but one precess- 
ing and one without any spin. At this mass ratio, both non- 
precessing and precessing binaries have observationally signif- 
icant power in higher modes; see, e.g.. Table V in Hannam 
et al. 

modes, with \m\ < I suppressed. By contrast, precession 
ensures that all five I — 2 modes are strongly excited, for 
all masses. 

For this simulation and reference mass (IOOM0) in par- 
ticular, antipodal decompositions always lose SNR for a 
significant fraction of orientations. One way to see this 
is by example, decomposing the waveform along z, J, Y 
in an antipodal decomposition, using each of the other 
combinations basis. None of these directions pro- 
vides an adequate antipodal basis for reconstructing the 
waveform along the others. More abstractly, antipodal 
decompositions fail when, even adopting a well-chosen 
frame (primed, below) to define the mode basis, a sig- 
nificant fraction of the orientation-averaged signal power 

m 

= (V'4;,„, V'4i„J = Plr-^ 

Im Im 

comes from harmonics other than the (2, 2)' and (2, —2)' 
modes. For this simulation and reference mass, the 
orientation-averaged SNR p for a source at 100 Mpc is 
~ 25. Using a (primed) frame aligned with V, at least 
nine modes contribute significantly (> 0.5%, summing to 
~ 99%) to p^; see Table HI In other words, for this sim- 
ulation and reference mass, other harmonics invariably 
contribute a significant fraction of the signal power, for 
many significant lines of sight. 

We have chosen a relatively low reference mass 
(IOOM0), for which higher-order multiples I > 2 produce 
a relatively small fraction of the orientation-averaged 
SNR (p). At higher masses, higher-order multipoles will 
produce even larger fractions of the total power [i^ and 



can lead to greater breakdowns of antipodal approxima- 
tions. 

To summarize, many simplifications extensively em- 
ployed for nonprecessing binaries break down completely 
for generic mergers. In general one cannot usefully ap- 
proximate the line-of-sight gravitational waveform with 
just two basis functions. As illustrated Table U for this 
concrete example, the |m| = 2 subspace is often but gen- 
erally need not be the most significant feature. For a 
nonprecessing binary dominated by the (2,2) and (2,-2) 
modes, the real and imaginary parts of the waveform 
along the z axis are proportional to the real and imag- 
inary parts of e^'*° and therefore are related by a sim- 
ple phase shift: the signal {ip4{z,^o)) produced by a bi- 
nary with initial orbital phase $0 is a phase shift by 7r/4 
(— 1-04(5, $0 + ''r/4)) of the signal produced by a binary 
starting with a similarly shifted binary phase. By con- 
trast, for a generic binary, multiple harmonics contribute 
significantly; for nearly no line of sight does the simple 
phase-shift relationship described above hold to the ac- 
curacy needed for gravitational wave data analysis. One 
cannot maximize over the orbital reference phase $0 us- 
ing an antipodal projection to extract coefficients of the 
"cosine and sine chirps" . The relevant maximum over 
coaelscen ce p hase can be found only if all harmonics are 
included gg. 



D. Constrained projections onto all harmonics 

More broadly, the signal from processing binaries can 
be projected onto a higher-dimensional basis. As only 
four degrees of freedom are needed to specify the line of 
sight (two euler angles, polarization, and distance to the 
source), the use of -/V > 4/2 basis vectors requires N — 2 
complex (2(A^ — 2) real) constraints, to insure the recon- 
structed signal estimate is consistent with what could 
generate it. 

Broadly speaking, constrained maximization methods 
have been extensively explored in gravitational wave data 
analysis. Constrained maximization methods arc pri- 
marily employed in the context of interpolating a signal 
manifold using a basis [48l . l49l . One particularly rele- 
vant example is Pan et al. [43|, who describe a strategy 
for detecting generic precessing low- mass binaries that, 
for each binary, finds the best-possible line of sight from 
the binary by first projecting data onto 5 basis signals, 
then solving constraint equations. In a similar fashion, 
the single orientation of a particular simulation that best 
fits a candidate signal can be recovered by projecting 
onto a simulation's harmonics. As our example illus- 
trates above, at least five and often more basis functions 
are required per simulation, with the number increasing 
rapidly with mass as higher harmonics become more sig- 
nificant in band j40| . However, in principle no obstacle 
exists to constrained maximization over source orienta- 
tion in an arbitrary reference frame. 

In practice, we expect that constrained maximiza- 
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tion of source orientation should be simplest in a frame 
aligned with the contours of source signal to noise (p^) 
and therefore in a frame roughly aligned with V, the 
principal eigendirection of {|L(qL{,)^^.^ In this frame, 
the hermetian overlap matrix Cab should generally have 
the simplest possible form 



{A,B) 



(20) 



where A is shorthand for normalized basis functions of 
time, proportional to a subset of tpii^{t). As a concrete 
example, the eigenvectors of Cab using the / = 2 basis 
states and our fiducial binary are fairly close to the basis 
states themselves. The reconstruction process for ip4{n) 
involves C explicitly and C~^ implicitly: 

Mn) = CABPB{n)\B} (21) 
PA ^ {A,^JH) (22) 

The four free parameters in the signal enter through ex- 
pressions involving C~^. First, the parameters pA must 
reproduce the total signal p^ = {ipA, ipi) along this line of 
sight, or 

2 



P 



Pa^abPb 



(23) 



The remaining three parameters correspond to specify- 
ing the orientation of the source, with three Euler an- 
gles. Two of the Euler angles can be specified using the 
principal axes of (LjaL^))^^, which has the form 

C 



{L(aLh)) 



M 



-1 
■'AB 

P" 



dflipA{n)* CaCbPB{n)) 



(24) 



The remaining Euler angle, the polarization angle, is con- 
nected to the phase of p. Though we adopt a multipolar 
basis in angle to define the basis set in time, the time 
domain basis functions \A) are not generally orthogonal, 
so pAin) are generally not proportional to spin- weighted 
harmonics.^ Conversely, only parameters in the corre- 
sponding 4-dimensional submanifold of the 2N (real) di- 
mensional space are consistent reconstructions of ipi from 
this source. 

Therefore, the overlap matrix C of the time-domain 
basis coefficients enters naturally into all reconstructions 
and constraints. Their underlying algebra will be sim- 
plest if this N X N complex matrix takes block diagonal 
form. However, our arguments above suggest that no 
fixed reference frame, not even one aligned with J, will 
make C block diagonal at all masses. 



^ For the purposes of this paper, we discuss constrained maximiza- 
tion of source orientation, assuming all harmonics arc known 
from simulations or post-Newtonian calculations. More broadly, 
by suitable interpolation across all masses and spins, a similar 
process of constrained maximization may allow waveform cata- 
logs to reconstruct all the basis waveforms of precessing binaries. 

^ In special cases such as nonprecesing inspiralling binaries, the 
time-domain basis functions are orthogonal and are propor- 
tional to spin-weighted harmonics. In this special case, the inte- 
gral on the right can be evaluated using representation theory of 
SU(2) H. 



IV. CONCLUSIONS 

Merging black hole binaries produce complicated 
orientation-dependent emission. For nonpreccssing bi- 
naries, the orbital angular momentum defines a pre- 
ferred orientation. For these simple binaries, a generic 
signal can be well-approximated by a linear combina- 
tion of 'ip4{L,) and '04 (—-£), or just one of the two func- 
tions tpii^L) away from the orbital plane. We demon- 
strate that these simplifying approximations break down 
for generic precessing binaries, by comparing different 
short, high-mass merger signals extracted along differ- 
ent orientations, using an extremely loose, data-analysis- 
motivated discriminator. Notably, our example shows 
that the final black hole spin angular momentum axis is 
emphatically not the natural emission orientation at late 
times. 

For nonprecessing and some precessing binaries, most 
detection- weighted gravitational wave signal power (p^) 
is radiated along the total angular momentum axis. Our 
counterexample also demonstrates that neither J, L, nor 
S generally corresponds to the direction of largest de- 
tectable gravitational wave signal. By contrast, the peak 
of SNR (p^) versus orientation is nearly along the domi- 
nant eigendirection of (^Lf^aLb)) ^ [Hi- 

For nonprecessing binaries, the orbit and kinematics 
define a preferred symmetry axis, relative to which the 
modulations encoded in the waveform appear particu- 
larly simple. For generic precessing binaries, particu- 
larly during the merger phase, different physical con- 
cerns (kicks; precession) lead to different preferred ori- 
entations. 

For data analysis purposes, we recommend an orienta- 
tion along the peak of the network-weighted SNR along 
the line of sight [p^) or the nearly- identical location iden- 
tified by the principal eigenvalue of (^L(^aLb)} j^^- This 
frame will change with mass. However, we anticipate 
PTF-style constrained maximizations over source orien- 
tations will generally be simpler in this frame. We will 
discuss the differences between this choice and the stan- 
dard frame in a subsequent publication. 

The breakdown of J as a preferred orientation comes 
as no surprise. As an extreme example, for sufficiently 
extreme-mass ratio binaries, the total angular momen- 
tum is always dominated by the black hole spin, carrying 
no information whatsoever about the motion. In fact, 
the relevant preferred orientation depends on how the 
detector "averages" the orbit. As a less extreme exam- 
ple, for BH-NS binaries, the orbital angular momentum 
direction is the preferred direction on any sufficiently 
short timcscale. For binaries whose masses and spins 
allow significant evolution in band, the relevant orienta- 
tion will become (roughly) a time-average of that vector 
[2H . Only if the orbital angular momentum precesses 
frequently about the total angular momentum in band 
will the total angular momentum arise as a preferred ori- 
entation. Even for high-mass-ratio binaries, precession 
can be irrelevant if it occurs out of band (e.g., before the 
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detected merger signal). 

Finally, in this paper we introduced several con- 
cepts useful for interpreting the signal from precessing, 
unequal-mass binaries. First and foremost, we describe 
a complex overlap and its natural maximization in time 
and polarization, a tool which efficiently eliminates two 
extrinsic parameters of generic signals even in the pres- 
ence of higher harmonics. Second, we emphasize that 
■04, being spin- weighted, naturally divides into two chi- 
ral parts; this decomposition corresponds to positive and 
negative frequency in the fourier domain. 



Appendix A: Chiral- weighted Orientation 

In general, the two extrema of need not be antipo- 
dal. Some methods for determining "the" optimal emis- 
sion direction treat n and — n symmetrically. To allow 
for different orientations of these two poles, we general- 
ize the tensor (^L(^aLb)} introduced previously to treat 
contributions from + and — helicities separately, if de- 
sired. For example, a matrix that sums up contributions 
from only right-handed emission is 



{L(aLb)) 



J^dfJcm4,i{f)L^aLb)Mf)/s^ 



rdfjdn\MfW/s, 



^1) 



This construct correctly identifies the asymmetric ex- 
trema of p± seen in even nonprecessing binaries; see, e.g., 
the red dot in the bottom panel of Figure [2j 



where we have required the two antipodal terms to be 
complex conjugates, as for a nonprecessing signal. In 
terms of these coefficients, we find the overlap P{n,fi') 
(before maximizing over time and phase) can be well ap- 
proximated by 



P 



22 I 



[Y*Y^ + Y1^YL2] 



O((-0422>'0422)) 

e^''^Y2{0)Y2{0') 



(Blc) 



-2i. 



^Y_2i0)Y_2i9') 



where we use the shorthand Y„i = Y^^^ (h) and simi- 
larly for y/„ to reduce superfluous subscripts and where 
in the last line we factor out the common e™*^ from 



^irn ■ above expression we assume rapidly oscil- 

lating, opposite-chirality terms like ('0422' V'422) cancel.^ 
We also make the trivial observation that the normal- 
ized overlap of a signal with itself is unity (P22 = !)■ 
This purely geometrical expression is largest when the 
two vectors are coaligned with their respective reference 
frames (i.e., when = 9' = 0). As an example, the spe- 
cial case of 6'' = has the form 



\P{z,n)\ 



Ds" e/2 



cos" e'/2 + sin^ e/2 



(B2) 



Appendix B: Complex overlap and antipodal 
decomposition for symmetric quadrupole-dominated 
emission 

When one or a handful of modes dominates, the an- 
gular dependence of the overlap P{n, fi') can be approx- 
imated everywhere via a handful of global coefficients. 
As a practical example, if 04 is dominated by symmetric 
I = |m| = 2 emission in the reference frame, the Wcyl 
scalar can be expanded "04 as 

04 = 0422W>"22"''(^)+V'422W*>"2"2'(^) (Bla) 
P22 = (0422,V'422)/|V'422l' = l (Bib) 
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Erratum : Is J enough? Comparison of gravitational waves emitted along the total 
angular momentum direction with other preferred orientations 
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The algorithm used to reconstruct tjj4.im to ^a{9, (f) in 
post-processing used a different convention for the spin- 
weighted spherical harmonics than the convention used 
to do the decomposition. Physically, this difference cor- 
responds to rotating the simulation results by 180° about 
the z axis. In this paper we performed all calculations 
starting with that codes' multipolar output. As a re- 
sult, all tensors and preferred directions derived from 
the original code were rotated by 180° relative to the 
simulation's frame. By contrast, the total initial and 
final angular momentum was correctly calculated, lead- 
ing to a discrepancy between the relative orientation of 
our preferred direction and the quantity we compared to. 
Figures 1, 2, and 4 must be replaced, as the relative ori- 



entation of J, L (simulation frame) versus the contours of 
^ and V (waveform frame) change. Figure 3, as a quali- 
tative illustration of the differences between the real and 
complex overlap, undergoes a trivial relabelling of the 4> 
axis. Table 1 remains unchanged. 

Figures 1 and 2 now show J is much closer to two key 
features extracted from outgoing radiation: first, our pre- 
ferred direction (blue line), which now precesses around 
J; and second, the largest values of |-04| and dE/dil. In 
Figure 4, the blue line showing the cumulative distribu- 
tion of P( J, n) now more closely resembles the black line 
showing the cumulative distribution of P{V*,n). 

The conclusions of the paper remain unchanged. 
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FIG. 1: Time-domain contours of \ip4,\ (solid) and dE/dtd^l 
(dashed), evaluated at the peak of the {l,m) = (2,2) mode. 
Both quantities are normalized to 1 at 6 = <^ = 0. Solid 
black dots indicate the location of the initial and final total 
angular momenta. The solid blue dot indicates the preferred 
axis identified by (1/(^1/;,))^ at this time. The dark blue path 
indicates how this preferred orientation evolves with time, 
several hundred M before and 30A'/ after this reference time. 
The top panel corresponds to the fiducial spinning binary. 
The bottom panel shows a binary with comparable mass ratio 
(mi/m2 = 4) but without any spin. 
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FIG. 2: Detector-weighted contours of p (dotted) and p+ 
(solid), evaluated at M — lOOi\/0, plotted versus the simu- 
lation's asymptotic spherical polar coordinates (i.e., aligned 
with the initial orbital angular momentum). Solid black dots 
indicate the location of the initial and final total angular mo- 
menta. The solid red dot indicates the preferred axis identi- 
fied by (^L(^aLb)) _^ [Eq. (??)] at this time; nearly identical 
preferred axes are identifies by (Lj^L;,))^^ and (L(„I/i,))^j._. 
As in Fig. [TJ the solid blue path indicates the evolution of 
(L(aib)) versus time. The dashed green path indicates the 
orientation of the (Newtonian, coordinate) orbital angular 
momentum oc r x v. The top panel corresponds to the fiducial 
spinning binary, while the bottom panel shows an nonspinning 
binary with comparable mass ratio {q = 4). In the bottom 
panel, the preferred orientation identified by {L(aLb)) ^ is 
slightly above the 6 = line; the two directions identified by 
i^L (^aLh)) ^ are not precisely antipodal. 
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FIG. 3: Area of high overlap: Fraction p of tiie spiiere 
for wliicii tiie overlap Pmax{v,fi) is above a given threshold 
P, versus the threshold, for v = z (blue), J (green), and the 
preferred axis of (^L(^a.Lf,)) (black). For comparison, the thin 
red curve shows the corresponding result for a binary radiat- 
ing only conjugate and orthogonal (2, 2) and (2, —2) modes 
and a reference axis v = z; see the appendix for a derivation. 
For this idealized nonprecessing binary, more than a third of 
all orientations have better than P = 0.95 (complex) match 
with the (2, 2) mode; all are well-fit with an antipodal decom- 
position. For the three reference orientations used here, the 
gravitational wave signal in that direction closely resembles 
the gravitational wave signal for only a small neighborhood 
surrounding that point. Results for simulations are estimated 
using ~ 10'' candidate viewing directions, each of which in- 
cludes alH < 4 contributions to ?/)4(?i). 



